arXiv:1502.01258vl [math.AP] 4 Feb 2015 


Enstrophy Cascade in Physical Scales 
for the 3D Navier-Stokes Equations 

Keith Leitmeyer 
February 5, 2015 


Abstract 

An enstrophy cascade is exhibited for the Navier-Stokes equations in physical scales 
independently of boundary conditions under physically reasonable assumptions on the 
flow. 


1 Introduction 

The Navier-Stokes equations Ut + {u ■ V)u = Aw — Vp, div u = 0, where u is the velocity 
vector and p is the scalar pressure, model the dynamics of an incompressible, viscous fluid. 
However, the physical theory of turbulence is based mostly on experimental observations 
and heuristic arguments. Two of the main features of Kolmogorov’s K41 phenomenology 
are existence of the energy cascade and scale-locality of the energy flux in 3D. The hrst 
rigorous mathematical result in favor of existence of the energy cascade was obtained in 
[9] in the setting of “stationary statistical solutions”, and the cascade took place in the 
wavenumbers. A related work [8] presented the proofs of existence of the enstrophy cascade 
and inverse energy cascade in 2D within the same mathematical framework (see also [14]). 
This approach concerns statistically steady-state turbulence, feeding off a non-trivial external 
force. Mathematical arguments conhrming the scale-locality of the flux were presented in 
[13, 6, 7], and [1] which features a rigorous proof of quasi-locality in the Littlewood-Paley 
setting. 

At least since G.I. Taylor’s fundamental 1937 paper. Production and dissipation of vor- 
ticity in a turbulent fluid ([16]), vortex stretching has been considered as a primary physical 
mechanism for creation of small scales in turbulent flows (cf. [10] for an overview of current 
efforts to establish rigorous mathematical framework for the vortex-stretching-anisotropic 
diffusion narrative). The tendency of a turbulent flow to self-organize in coherent vortex 
structures, most notably in vortex hlaments, has been well-documented in direct numerical 
simulations, see, e.g., [15, 17, 12]. Strong local anisotropy and quasi low-dimensionality of 
turbulent flows in the vorticity description points to plausibility of existence of 3D enstrophy 
cascade. However, in order to efficiently exploit depletion of the nonlinearity caused by the 
local anisotropy, and in particular, by local coherence of the vorticity direction (cf. [3]), it is 
necessary to be able to formulate dynamics of the turbulent cascades directly in the physical 
scales of the flow. 
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A spatial multi-scale ensemble-averaging process designed to detect significant sign- 
fiuctuations of a density in view, at a given physical scale, was introduced in [4], This 
was then utilized, taking the density to be the energy flux-density, to obtain a proof of 
existence of energy cascade and scale-locality of the flux in physical scales of 3D turbulent 
flows, in the case of free (decaying) turbulence. A subsequent work [5] presented a proof 
of existence of 3D enstrophy cascade in a cylinder in space-time under the assumptions on 
the geometry of the flow and smallness of a Kraichnan-type micro-scale, as well as several 
technical assumptions. One technical assumption, that the enstrophy - localized to the 
macro-scale domain of interest - was smaller than a given constant, was (as pointed by the 
authors) less than satisfactory, both from mathematical and physical viewpoints. 

In what follows, we make modifications to the proof in order to replace this assumption 
with a much weaker assumption that can be thought of as a bit more restrictive version of the 
Kraichnan-type condition, and is consistent with expected spatial complexity and uniform 
locality of fully developed turbulent flows. 

The key new ingredients are redesigning the ensembles of test functions in the construc¬ 
tion of the ensemble averages in order to establish an explicit relation between the ensemble 
averages at two different scales, and formulating uniform locality in terms of Morrey-type 
quantities. 

The main theorem states that the ensemble average of enstrophy fluxes within a range 
of scales is comparable to the modified mean palinstrophy (in particular it is positive). In 
the first section the ensemble framework is presented. Next, we state all of the assumptions 
on the flow used to prove existence of the enstrophy cascade. Then the proof of the main 
result is given. 

2 Ensemble Average Framework 

The inward enstrophy flux through a sphere is given by — ■ nda, where n is the 

outward unit normal vector to the surface and the enstrophy is the norm of the vorticity 
u = curl u. By Stokes’ theorem, — ■ nda = f^(u ■ V)lj ■ u dx. Instead of the 

sharp cutoff at the boundary of 5, we need something smooth. We will be working with 
the quantities f (u ■ V)lj ■ tfju dx = — f dx, where is a smooth function equal 

to 1 on 5, supported on 2B, and with inward pointing gradient, with certain bounds on its 
derivatives. 

Definition 1 (Refined test function). A {Co, p) test function at scale R is any tjj G C°°(M^) 
supported in a ball of radius 2R with 0 < ■^ < 1, |V'^| < , and |A-^| < 

Definition 2 (Ensemble). (Parameters R,Co, p, Ki, K 2 ,'ipo,Ro) a test function fjQ = 
'ipo,Ro- The region of interest is B{0, 2Ro). 

An ensemble is a collection of {Co, p) scale R test functions satisfying the following 

properties: 

1 . fji ^fjo 

2. {Ro/Rf <n< Ki{Ro/Rf 
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3. No point of B{2 Rq, 0) is contained in more than K 2 of the supports of ipi. 

Choosing larger Ki, K 2 allows ensembles with greater global and local multiplicity re¬ 
spectively. 

Definition 3 (Ensemble average). For a function f, denote by {F)r the ensemble average 
^ YTi=i w I dx, and Fq = ^ f ffjo dx. 

Property 1 above is needed to compare (E)/j to Fq. Note that test functions near the 
boundary of the support of fjQ will have small integrals (due to Property 1), effectively 
skewing the ensemble average towards zero. Larger Ki, K 2 allow ensembles that have higher 
weight on functions away from the boundary, making the skewing insignificant. 

These ensemble averages can be viewed as a way to detect whether a function is signifi¬ 
cantly negative at some spatial scale. If every ensemble average {F)^ (for fixed parameters) 
is positive, no matter how one arranges and stacks the test functions, then the function is not 
significantly negative at scale R. Increasing Ki and K 2 lowers the threshold for a function 
to be considered significantly negative. 

Many ensembles can be constructed by applying Lemma 2 to ipQ and varying the mul¬ 
tiplicity of the resulting functions (Assume -ipo satisfies the stronger CQ-bounds to get an 
ensemble with Co-bounds). 

The following lemma states that ensemble averages (at any scale) of positive functions 
are comparable to the large scale mean. The proof immediately follows from the definitions. 

Lemma 1. /// > 0 then -^Fq < {F)ii < K 2 F 0 . For slightly modified ensemble averages, 
we have lYJi w ! dx <-^ 2 -^ f /^o dx (d > 0). □ 

Using a refined partition of unity, one can turn larger scale ensembles into smaller scale 
ensembles. 

Lemma 2. Any (Co,p) scale R test funetion is a sum 0 / 64(i?/i?')^ {Cq, p) scale R' test 
funetions (where C( depends only on Cq, R > R!). 

Therefore for all {Ki, K 2 ,Co)-ensembles at scale R and every R' < R, there exists a 
(647^1, 87^2, Co)-ensem6/e at scale R' such that {F)ii = {F)ri. 

Proof. Let be a (Cq, p) scale R test function. Now to construct the partition of unity, take 
a (Co,p), scale R' 3D test function qq, centered at zero and equal to 1 on [—R,R\^ (such a 
function exists as long as Co isn’t too small). Define Pp = go{x — 2R'p), where p E If. Then 
1 < Ep S'? < 2 so we may define hp = Pp/ Y.q 9q- 

Some calculus shows that |Vhp| < ^§rhf and \Ahp\ < so |V('^hp)| < 

’^{fhpY and \N.{fjhp)\ < ('^hp)^^~h Fewer than 8\R/R'Y < 64{R/R'Y of the 

functions fhp are nonzero, and for any x, 'fp{x) Y 0 for at most 8 functions. 

Since ip = Ept^^P) have the first claim. For the second claim, given an ensemble 
{ipi}i, the new ensemble will be {ipihp}i^p. □ 
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3 Enstrophy Cascade 

3.1 Assumptions 

Let be a domain in and u a Leray solution to the Navier-Stokes equations on hi x (0, cxd). 
The enstrophy cascade will occur on a cylinder 5(0, 2i?o) x (0, T) where 5(0, C hi, 

5o < 1 and T > Rq {Rq will be the macro-scale in the problem). Centering 5(0, 25o) and 
starting the time interval at zero is for notational convenience. 

It is required that the solution u has supjg(o,r) Ib{o2Ro+r^^^) dx < oo. If hi = R^ 

we require supjg^Q j.) \uj{xR) \ dx < oo. This is guaranteed if u is a classical solution on 
(0,T) with IIcuoIIli < oo, or is a Leray solution with finite Radon measure initial vorticity 
given by the retarded mollification method used in [ 2 ]. 

Fix a {C,p) scale 5o test function 'ipQ for spatial localization. Let a temporal cutoff 
function rj G C°°[0,T] be such that \ri'\ < ^ 77 ^, 0 <? 7 < 1 , -/7 = 0on [0,T/3), and 77 = 1 on 
(2T/3,T], where T > Rq. For any test function ip, dehne <p{x,t) = ip{x)ri{t). 

In [11] it is shown that Assumption 1 implies that the localized enstrophy is bounded on 
(0,T). Using the localized Biot-Savart law, this implies a bound on \ \u(t)\\Hi(B{Ro))- Then by 
the partial regularity theory of the Navier Stokes equations u is smooth on (0, T] x 5(0, Rq). 

For simplicity we consider hi = R^ only. If hi is a bounded domain, the localized Biot- 
Savart law used in [ 11 ] has a number of lower order terms, which will require a small modih- 
cation of the assumptions. Namely, that in Assumption 2 and R = ((Jo//9)^/"^ 

in Theorem 2. 

Assumption 1. Let ^ = oj/\uj\ he the vorticity direction field. Assume there exist M,Ci 
such that I sin79(,^(a;,f),.^(7/,l))| < Cpx — for a.e. {x,y,t) in (5(0, 25o) H {|Vm| > 
M}) X 5(0, 25o + X (0,T), where (p{zi,Z 2 ) denotes the angle between the vectors Zi 
and Z 2 . 

This assumption is based on numerical simulations which suggest that regions of intense 
vorticity self organize into coherent vortex structures and in particular, vortex hlaments. 

Assumption 2. Denote the scale-Ro mean enstrophy by Eq = A J ^ J h^uj\‘^(j)Qp~^dx dt, the 
modified mean palinstrophy by Pq = A J ^ J |dxdl + / \^{x,T)\‘^fi’odx, and the 

modified Kraichnan scale by ag = R is required that ao < fiRo (t) < (5 < 1 is a 

constant depending only on Cq,Ci, M, Ki, K 2 , and Bt := supjg^Q j’) ||ca||Li(n), and fd shrinks 
to zero as any of them increase to infinity). 

The second term of the modihed palinstrophy arises from the shape of the temporal cutoff 
77. 

Denote by = M^’'^(5(0, 25o + -^o ^)) Morrey space of functions / such that 
^^^y,R h(y,R)nBiopR,+Ryfi ^ ^ ' 

Assumption 3. Assume u(t,x) G L^(0,T; with where 

C depends only on fd,Co, Ki, K 2 . 
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Assumption 3 will be used with the bound | |i^| |L2((o,T)xB(a;i,2i?+R2/3)) < |a;| |^ 2 ^ 2 , 5 . 

All Leray solutions have ||ci;|| 2 , 2£,2 < oo. We need slightly more - that | |a;| |^ 2 jy^ 2.9 < oo. 
Assumptions 2 and 3 will be true if ||a;|L 2 ^ 2,5 is sufficiently small relative to Pq, that is, we 
consider high (time-averaged) spatial complexity of u in B{0,2Rq). 

The hnal assumption used to prove the enstrophy cascade is that the enstrophy doesn’t 
drop off too much at time T. 

Assumption 4 (Modulation). J \u{x,T)\'^^jJo{x) dx > J \u{x,t)\'^ipo{x) dx. 


3.2 Theorems 

To work with the enstrophy we use the vorticity form of the Navier-Stokes equations: 


ut + {u ■ V)u = {u ■ V)u -|- Au. 


To use the ensemble average framework, we take as our function / = — T jJ'(u-V)u-u rj dt 
so that J fipi dx = ^ f ||a;p(u ■ V0i) dx dt. This time averaged enstrophy flux along 
will represent the amount of enstrophy flowing into the scale R if 0^ is taken to be constant 
on B{xi,R) and with inward pointing gradient. 

The following lemmas will be used in the proof of the enstrophy cascade, and are essen¬ 
tially as they appear in [4, 5]. 

Lemma 3. For a solution to the Navier-Stokes equations u that is smooth on [0,T] x 
B{xi,2R), 



1 

2 


u:\^{u ■ V(j)i) dx dt 



u{x, T)\‘^^jJi{x) dx -f 



I Va;|^0i dx dt 



-\u\^{dt(t)i -f A0j) dxdt 



j{u ■ V)u ■ (piU dx dt. 


Proof. Using integration by parts and that u is divergence free, ||a;p(M-V0i) = —{u-'V)u-(f)iU. 
From the Navier-Stokes equations, — (u ■ V)a; = dtu — Au — {u ■ V)u. Integrating in space 
and time against (piU yields — fj' f (u ■ V)u ■ piU dx dt = f dtu ■ piU dx dt — J Au ■ 
piU dx dt — fJ' f (u ■ V)u ■ piU dx dt. Now 


dtU ■ PiU dx dt = 


1 1 

-^dt{u ■ piu) - -u ■ {dtpi)udxdt 


= / \u{T)\^pidx - 


\u\'^dtpi dx dt. 


and 


Au ■ pjU dx dt = 


|Va;|^0i dxdt — 


-^\u\^Api dx dt, 


simply by integration by parts and the fundamental theorem of calculus. Putting all of the 
equations together completes the proof. 

□ 
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Lemma 4. For a divergence free function u G with u := V x u, we have {u ■ V)u ■ 

cj(x) = cF.V.fcj(x) X cj(y) ■ G^(x,y)dy for a.e. x, where iG^{x,y))k = 

= + j^Ukix). 

Proof. For any divergence free Schwartz function u, we have Au = —curl ca, so n = c |i| * curl 
u, and diU^ = cdi^ * Cjkidki^^ = cdiduy^ * e^kii^K Then 


uj^djU^uj\x) = coj^oj\didk -.—r * ejkiuj^){x) 


cP.V. 


j Cj(x) X Cj(y) ■ Gu,(x, y) dy. 


By density of Schwartz functions in the following holds in L^: 


{u ■ V)n • u = cP.V. 


u{x) X u{y) ■ Guj{x,y) dy. 


□ 

Now we wish to bound the vortex stretching term by integrals of positive functions in 
order to use Lemma 1. 


Lemma 5 

f 

'0 


{uj ■V)u- (fiudxdt < c||a;||*^sup / -^\u:{x.,t)\^ipi{x) dx + / |dxdfj + 


c' + c"||a;| 
^2 


\u\^(j)l'^ ^dxdt 


(I leal I* := ||a;||i 2 (B( 2 ;,,ij)x(o,T)), constants depend on M,Bt,Go). 

Proof. Jp J{u ■ V)u ■ (fiUdxdt = /o^/{|v„|<M}(^ ■ • 4>i0Jdxdt + ■ 

(piudxdt. The hrst term is easily bounded by ^ J |cap0^^~\ For the second term, the 
Biot-Savart law gives us 


'0 a{|v«|>M} 

r/ 

'0 a{|Vn|>M} 

f! 

'0 a{|Vn|>M} 


(ca ■ V)n • (fiOJ dx dt = 


P.V. / ca(a;) x u){y) ■ G^^^x, y)(l)i{x) dy dx dt = 


ao J{\^u\>M} J{\x-y\>B?/^} 

The second term (4) is bounded by 


P.V. / ta(a;) X uj{y) ■ Gi^{x, y)(i)i{x) dy dx dt 

'{|x-j/|<_R2/3} 

ta(a;) x u{y) ■ Guj{x,y)(f)i{x) dy dx dt. 


1 


( 1 ) 

( 2 ) 

(3) 

(4) 



{\x-y\>R'2/3} 


ca(a;)|^|ta(|/)|0da;) d|/dxdt < — sup ||ca(t)||i;,i / / |ca|^0^^ ^dxdt. 
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For the first term (3), since 

|a;(a;) x u{y) ■ G^{x,y)\ < |a;(a;)||a;(l/)|| sin(^(^(a;),^(l/) )||G^(a;,|/)| < 


we have 


'0 J{\Wu\>M} 

pT r 


/ u{x) X uj{y) ■ Gi^{x, y) dy 4)i{x) dx dt < 

'{|x-j/|<_R2/3} 

\u{y)\\u{x)W,. 


'0 J{\Vu\>M} J{\x-y\<R'^/^} \x — 


-(j)i{x) dy dx dt < 


C / ll'^llL2(B(a;i,2_R+_R2/3)ll 10 / '^ni3/2'^^ < 

Jo 

pT 

c / ll^^llL2(B(x,,2ij+ij2/3))||0y^x;||2||V(0j/^a;)||2rft < 

Jo 

1 j‘2i 

c||a;||,(^sup||^y^a;||^ + ^ \\V{(j)y^uj)\\ldt^ < 

c||i^||*(^^sup ||'?/’y^a ;||2 + 2 f f \Vu\‘^(pi dxdt + — 


\u\‘^(pf^ ^dxdt^, (11) 


using |V(0y^a;)p < 2|Va;p0j + [cap < 2|Va;p0i + for fhe last inequality 

(11), and Hardy-Littlewood-Sobolev to reach (7), Gagliardo-Nirenberg to reach (8), and 
Cauchy-Schwarz to reach (9). Collecting the bounds on the various terms proves the lemma. 


Theorem 1. j^Po < {P)r < {K2 + j^)Po for R = Cq/IJ, and {Ki, K 2 )-ensemble averages. 
Proof. For an individual test function we have 


P. : = 


^ 1 /* 1 P 

j -^\uj\^{u ■V(j)i) dx dt = ^sup j -^\uj{x,t)\^ipi{x) dx + j j \Vu:\^(j)idxdt 

|a;p((9t0i -|- A0j) dxdt — f [{uj ■ V)m • (piudxdt =: Ai — Bi — Gi. 


Using Assumption 4 and Lemma 1, i^Po < nTi?3 YfJl — P 2 Po- Next, \Bi\ < 

W lo I dx dt so < ^-^0 < ciF 20 ^Po < ^Po for an appropriate 

choice of /3. 

Using the vortex stretching term lemma and Assumption 3, 

1 ^ PC 

\-iij^EA<(‘= + JM)^Eo + c"\MK2P» < 




^llL2M2’'jf^2Po < ^^^0 
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Then {F)r = A - Bi - Ci hence ^Pq < 

{F)r < {K, + ^)Po 

□ 

Theorem 2 . 25 ^< {F)r < { 8 K 2 + 2 ^)^o for {Cq, p, Ki, K 2 )-ensembles with a^ffi < 
R ^ Rq . 

Proof. By Lemma 2, every (Cq, p, Pi, P 2 ) scale R ensemble average is equal to some (Cq, p, 
64 Pi, 8P2) scale ao//9 ensemble average, which satishes the desired inequalities. □ 

In particular this holds for ensembles of test functions with inward pointing gradient that 
are constant on B{xi, R). Thus in this precise sense solutions of the Navier-Stokes equations 
that satisfy the assumptions exhibit an enstrophy cascade. 


4 Conclusion 


Fix Rq, T, Ki, K2, Cq, Ci, M, Br, and q. There will be an enstrophy cascade on P(0, 2Po) x 
(0,T) for any Leray solution with Eq/Pq and \ \oj\\T2j^^,q / Pq sufficiently small, as long as there 
is coherence of vorticity direction where |Vm| > M and the macro-scale enstrophy before 
time T never exceeds twice the local enstrophy at time T. The assumptions are consistent 
with strong local anisotropy, spatial complexity and uniform locality of fully developed 3D 
turbulent flows. 
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